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In this paper we revisit the onset of the instability of the solid state in classical systems within 
self-consistent phonon theory (SCPT). Spanning the whole phase diagram versus volume and versus 
pressure, we identify two different kinds of mechanism: one mainly relevant at constant volume, 
associated with the vanishing of the SCPT solution; and one related to the disappearing at a 
spinodal temperature of the solid phase as a metastable energy minimum. We show how the first 
mechanism occurs at extremely high temperatures and it is not reflected in any singular behavior 
of the thermodynamical properties. In contrast, the second one appears at physical temperatures 
which correlate well with the melting temperature and it is signalized by the divergence of the 
thermal compressibility as well as of the the lattice expansion coefficient. 

PACS numbers: 64.70.dm, 63.10.+a, 68.35.Rh, 62.50.-p 



I. INTRODUCTION 

The solid-liquid transition is one of the most common 
and extensively studied phase transition in condensed 
matter. As a first order transition, overheating the solid 
beyond the melting temperature is possible in nature,—"— 
following the metastable equilibrium state (Fig. [T]). 

This route is however limited by the intrinsic instabil- 
ity of the metastable phase one is considering. Focus- 
ing on the metastable solid phase, a priori the maximum 
temperature of overheating (for a surface- free and perfect 
crystal without any defect or impurity) could be identi- 
fied with the temperature above which the solid phase 
is not a sustainable phase. The breakdown of the solid 
phase independently of the competition with the liquid 
phase is thus associated with the concept of intrinsic in- 
stability of the solid phase, i.e. with the breakdown of 
the conditions which make a solid phase sustainable even 
at a metastable level. 

Historically, the problem of the mechanical melting^ of 
a solid was discussed indirectly by the one-phase theories 
addressing the standard meltingii^rJ^ The most two fa- 
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Figure 1: (color online) Pictorial evolution of solid/liquid 
Gibbs free energy G as a function of temperature T . At the 
melting temperature Tm , the Gibbs free energies of the solid 
and liquid phases are equal, the system undergoes a first or- 
der transition. A metastable solid phase (dashed line) can be 
still obtained up to a maximal point Tmax. 



mous theories are the well-known Lindemann criterion^ 
and the Born rigidity catastrophe Although these the- 
ories cannot explain the standard melting, they have been 
reconsidered, with a modern point of view, for a theoret- 
ical description of the instability of an overheated solid. 
Such a phenomenon, which was once thought to be un- 
obtainable, has become now practical due novel experi- 
mental advances in heating techniques and in the fabrica- 
tion of special samples, as it has been reported in a wide 
number of different systemsii"— These experiments have 
stimulated a renewed theoretical interest for the problem 
of the mechanical melting. 

Recent molecular dynamics simulations have been also 
carried out to clarify the underlying microscopic mech- 
anism which sets the stability limits of a overheated 
solid.—"— In particular, Jin and coworkers found that 
above the equilibrium melting point T™, local lattice in- 
stabilities governed by both Born and Lindemann criteria 
occurs at a well defined mechanical melting temperature 
7max,~ at which the material cannot survive in crys- 
talline order for any finite time interval. Interestingly, 
the temperature of the mechanical melting was also as- 
sociated with a sudden and drastic rise of the atomic vol- 
ume corresponding to a large peak of the compressibility 
of the system.— 

The crystal instability comes ultimately from anhar- 
monic effects that may soften the lattice when the ther- 
mal fluctuations are large. To this issue, the mechani- 
cal melting was also addressed in the past by using ex- 
tensively the self consistent phonon theory (SCPT),— "— 
which represents a suitable method to take into ac- 
count the anharmonic effects of the atomic oscillations. 
Within this microscopic approach, the breaking of the 
self-consistent solution of the SCPT is thought to repre- 
sent the mechanical instability of the solid phase and it 
defines a maximum temperature T* which was commonly 
assumed as upper estimation of the melting temperature 
(the variational free energy obtained by the SCPT is an 
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Figure 2: (colors online) Schematic phase diagram for the 
breakdown of the mechanical stability of the solid phase 
within the SCPT. Increasing the temperature at fixed vol- 
ume (solid green vertical line), the solid phase disappears at 
T* which represents, in this case, the highest temperature at 
which the variational approach (SCPT) has a self-consistent 
solution (red line) . Due to the lattice expansion, the mechan- 
ical instability at fixed (reasonable) pressure {P < Pc) is first 
encountered in terms of the disappearing of the metastable 
solid state, defining a spinodal temperature T — Ts{P) (blue 
long dashed line) within a well-defined SCPT solution. Below 
the line Ts{P) the compressibility kt is negative. 



upper bound of the exact free energy). However the val- 
ues of the temperatures T* obtained in this way are much 
higher than the experimental melting temperatures 
(even of two orders of magnitude) observed in several 
systems 1^ Moreover, they compare badly also with the 
experimental overheating temperatures Tmax which are 
generally located between T,„ and ^ l.5Tjn^ 

In this paper we revise the self-consistent phonon the- 
ory focusing on its extension for crystals at constant 
pressure P, taking into account in particular the lat- 
tice volume variation as a function of both temperature 
and pressure. Our main results are summarized in the 
schematic phase diagram depicted in Fig. [2l 

In particular we find that within the SCPT the phys- 
ical description of the mechanical instability of the solid 
phase approached at constant volume is qualitative and 
quantitative different from the mechanical instability 
achieved at constant pressure. In the first case, the insta- 
bility of the solid phase occurs at rather high tempera- 
tures and it is associated with a breakdown of the SCPT 
solution at T*, not reflected in any singular behavior of 
the thermodynamical quantities. On the other hand, the 
mechanical instability approached at fixed pressure can 
be properly interpreted in terms of a spinodal temper- 
ature Ts at which the solid phase, always defined, dis- 
appears as a metastable minimum. In the latter case 



the mechanical instability is thus reflected in the diver- 
gence of the thermal compressibility kt as well as in 
the divergence of the lattice expansion coefficient a at 
T = Ts- This means that small fluctuations of the two 
controlling fields (T, P) produce large volume fluctua- 
tions, AV/V = aAT and similarly AV/V = -ktAP, 
pointing out that the system is mechanically unstable, 
as previously discussed phenomenologicallyJ^ 

To show the main features of our analysis, we focus 
initially on the specific case of solid Argon and we gen- 
eralize later these results for other rare-gas solid systems 
(Ne, Kr, Xe), resulting in Ts m fair agreement with the 
experimental values of the melting temperatures T^. 

The paper is structured as follows. In Sec. |lT] we 
present our model for generic classical solids, recalling 
the SCPT for fixed volume and discussing its extension 
for the fixed pressure case. The results for the specific 
case of solid Argon are shown in Sec. Illli pointing out 
how it is possible to achieve two qualitative and quanti- 
tatively different mechanical instabilities by working at 
constant volume or at constant pressure. We also dis- 
cuss there the singular behavior of the thermodynamical 
properties of the system. In Sec. llVl we explain the phase 
diagram shown in Fig. [5J In the last Section|V]we drawn 
our conclusions. 



II. THE MODEL 

We consider a classical solid formed by N particles with 
mass m and interacting via an isotropic U (r) = U{r) pair 
potential: 

where i ^ 1, . . . ^ N and where Pi and are the momen- 
tums and the positions of the particles. The classical 
partition function of the system reads: 

ZiT,V)^JjYl^ e-«/^-^ (2) 

i 

which is a function of the total volume of the system V. 
The statistical average for a generic physical quantity O 
is given by: 

(O)-|//n^O(r.,p,0 e--A^-. (3) 

For sake of clarity, we consider here a monoatomic crys- 
tal with cubic symmetry. At zero temperature the crys- 
tal is frozen and, in absence of any defect and far from 
boundary surfaces, the particles are fixed at their lattice 
positions — with a ~ |Ri — Rj| when (i,j) are 
neighboring atoms. Because we work at fixed particle 
number, we have the simple relation V oc Na^. 
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At low temperature, the thermal fluctuations can be 
described in the harmonic approximation where we ex- 
pand the interatomic potential for small lattice displace- 
ment of the particles = — around their average 
position R,;. By assuming a mean field approximation 
(the Einstein model), we can write: 



where Rij is a short notation for Rij 



Ri — Ril and 



k is the elastic tensor responsible for the restoring force 
on each particle. In the cubic symmetry the tensor k is 
diagonal and isotropic in the three axial directions a = 
x,y,z, so that kap = kSap. We thus obtain the well- 
known result for the mean thermal fiuctuation (m^) = 
SksT/k = 'iksT /mujQ where the Einstein frequency is 
cjQ — \fkjm. As we raise the temperature to approach 
the melting point and beyond it (the overheated regime), 
anharmonic effects are expected to be relevant. We can 
take them into account by using a variational method, 
the self consistent phonon theory, which we discuss in 
the next paragraph. 



A. The SCPT at fixed volume 

We recall the standard scheme of the variational meth- 
ods. A (quadratic) trial Hamiltonian i?v containing vari- 
ational parameters is introduced. The corresponding par- 
tition function reads: 



zAT.y) = // n 



d r,dp, 

/l3 



(5) 



Then the Gibbs-Bogoliubov inequality gives as upper 
limit of Z: 

Z = Z,(c-(^-^^)/'^--^), < Z, c-««-^^)/'^--^>v . (6) 

where {■ ■ ■)^ denotes the average as in Eq. ([3]) on the trial 
Hamiltonian H^. From Eq. ([6]), we have the following 
inequality for the Helmholtz free energy F{T, V): 

F = -kgT In Z < -keTlnZ^ + {H - H^), = , (7) 

where the r.h.s is the variational free energy Fv{V,T). 
Then the variational parameters are determined by min- 
imization of Fv 

We now apply the above described variational ap- 
proach within the context of the Self Consistent Phonon 
Theory (SCPT) assuming a harmonic local model for any 
temperature. We write 



i/. = i^L/(R.,) + 5] 



2m 



(8) 



where the local force constant k^ is the variational param- 
eter. It describes an effective elastic force with frequency 



LUy = ^Jkyjm. We use now the hamiltonian , Eq. ([8]) , 
to calculate the r.h.s. in Eq. ([7]), the variational free 
energy Fv{T, V) whose the first term reads: 



1 / k T 

fcisTlnZv = -^J7(R.y) - STVfcsTln ' ^ 



, (9) 



By using the cquipartition theorem for the quadratic 
term, we have also directly the average of the exact 
Hamiltonian over the trial Hamiltonian: 



and also the average of the trial Hamiltonian: 
{H.)y^lYlU{Il,,) + 3NkBT. 



(10) 



(11) 



Summing up Eqs. and ([TT|) . we finally get the 

variational free energy per particle: 

^ = -3fcBTln ( M V— y {Uir^ ~ r,)) --fc^T . 
N \fujyl 2N^^ ^"^ 2^ 

(12J 

It is convenient to introduce the smeared potential U 
defined as: 

U{B..,,,ul) = (C/(r,-r,))^ 

= /_^[/(k)e'i^-Ve'>»-"^)'' 



where ur, 



i2nr 
drU{r + R,j) ■ 



SkeT 



-r7(4/3)< 



'4 9\3/2 



(13) 



(14) 



Eqs. P^ - (fTi|) define at this stage the explicit form of 
the variational free energy F^ which has to be minimized 
with respect to fcv. Assuming isotropy along the x, y, z 
directions, from the condition dF^/dk^ = we get: 



-E 



A^^./(27r)='" 3 



E E 



(15) 



which must be solved self consistently since the smeared 
potential U depends implicitly on k^ via Eqs. (fT3)) . (|14p . 
In the case of a fully isotropic potential f7(Ry ) = U{Rij), 
and reminding Eq. p^ . we get thus the compact self- 
consistent solution for u?- 



3kBT 




d^U{x,ul) ^ 2 dUjx, ul 
dx^ X dx 



x—Rij 

(16) 
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The SCPT allow us to evaluate in a self-consistent way 
the mean thermal fluctuations (and hence the ef- 
fective elastic constant fcv) for given temperature. In 
the literature such an approach has been employed 
in the analysis of lattice mechanical stability against 
anharmonic fluctuations as temperature raises at fixed 
volumci^ The temperature of instability was thus iden- 
tified with the maximum temperature T* above which a 
self-consistent solution for disappears. 



B. Extension at fixed pressure 

After the minimization, the SCPT provides us an 
approximated expression for the Hclmholtz free energy 
F(T,V) ~ Fv(r, as function of the volume V. In 
order to describe how the volume per particle V/N 
evolves with temperature and pressure, we introduce 
in the usual way the Gibbs free energy G defined as 
G(T,P) ^ F^[T,V{T,P)] + PV{T,P), where now the 
volume V — V{T, P) has to be considered itself as a 
function of temperature and pressure, and it is given by 
the relation [dG / dP)rp ~ V or, equivalently, by inverting 
the relation 



dF^{T,V) 



dV 



(17) 



which gives an implicit definition of the volume V 
ViT,P). 



III. RESULTS 

In order to keep the calculations at the most analytical 
level, we assume in the following a particularly suitable 
form for the potential U (r) which reproduces the stan- 
dard molecular potentials in the whole effective range for 
distance experienced by the particles at zero temperature 
as well as in the vicinity of the instability temperature 
Ts- We consider in particular a potential given by the 
linear combination of two Gaussian representing a long 
range attractive tail and a short-range repulsion)^ 



C/(r) 



Uo 



/?+e 



(18) 

where r = |r| and where Uq is the potential minimum at 
the point r ~ r^. By a proper choice of the two parame- 
ters /3+ , /?_ it is thus possible to reproduce the behavior 
of the standard molecular potentials as, for instance, the 
Lennard-Joncs or the Morse potential in the range which 
we are interested. We stress that the two parameters 
/3+,/3_ are independent on temperature. An example is 
given in Fig. [31 where we compare the Gaussian-like po- 
tential of Eq. ([rs]) with the Morse potential describing 
Argonj^ 

In Fig. [3] the lattice parameter a (black dots)and the 
range of the corresponding lattice fluctuations (double 
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Figure 3: Comparison between the Morse and the Gaussian- 
like potential for the case of Argon. Model parameters are 
here ro = 3.76 A, Uo = 146.8 K and 7 = 1.55 A"^ for the 
Morse potential, and ro = 3.76 A, Uo = 146.8, P+ = 0.45 
and j3- = 0.16 for the Gaussian-like model, obtained 
by the best fitting of the Morse potential. The filled circles 
and horizontal double arrows represent the lattice parameter 
and the corresponding lattice fluctuations at low temperature 
(T ~ 1 K) where a ~ 0.98ro and at the instability tempera- 
ture T ^ Ts ^ 227 K where, for P = 0, we have a ^ 1.2ro 
due to thermal expansion. 



arrows) are also shown for T w and for T ^ Ts{P ^ 
0) ~ 227 K as representative cases of low and high tem- 
perature respectively. As we see, the agreement between 
the Morse potential for Argon and our Gaussian-like po- 
tential is practically perfect in all the physical range of r 
spanned at both low and high temperature. This agree- 
ment is particularly important since it assures that the 
results of the variational method for two potentials are 
indistinguishable. 

For the Gaussian-like form of Eq. (|T8)) . the smeared 
potential reads: 



C/(r) = 



X exp 







13+ - 13- ^ 
r 



.=± (1 + paUl) 
2 



3/2 



1 



- r, 



(19) 



In the zero temperature limit (m^ = 0) the smeared po- 
tential reduces to the bare potential U ~ U and the elas- 
tic constant Eq. to the bare elastic constant = kg. 
The zero temperature/zero pressure nearest neighbors 
distance oq = a(T = 0, P = 0) is found by minimizing the 
classical ground state energy Eq = / "2) J^i^j U (Rij) . 
For the Argon parameters considered in this section we 
get fco ~ 350.5 C/o/7'g and oq ~ 0.98 7'o which is in good 
agreement with the experimental value Uq^^' = 0.99 tq.— 
Note that the value gq ~ 0.98 ro is slightly different from 
the potential minimum ro due to the contribution of the 
second nearest neighbors atoms in the total energy. By 
increasing the temperature, both the thermal fluctua- 
tions and the lattice parameter a will have a non trivial 
dependence on temperature. 
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A. Constant volume 

We discuss first the case where the lattice parameter a 
is assumed to be independent of T and set to the value 
a = oq, i.e. its value at zero temperature and zero pres- 
sure. In this case, the effective elastic constant fcv(T) and 
the lattice fluctuations u^iT) can be simply obtained as 
functions of the temperature by the self-consistent solu- 
tion of Eqs. p^ - ([TC| . setting the a = cq in the lattice 
sum. The graphical solution of Eq. (|16p for three repre- 
sentative temperatures is shown in Fig. m 




Figure 4: Graphical solution of the self-consistent equation 
Eq. (|16p for fixed lattice parameter a = ao and for the Argon 
parameters. The dashed lines correspond to the l.h.s. of 
Eq. (|16|) for three representative temperatures while the solid 
line represents the r.h.s. of Eq. (|16|l which does not depend 
parametrically on the temperature. 
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Figure 5: (a) Mean lattice fluctuations u1{T), (b) effective 
elastic constant kv{T), solutions of the self-consistent equa- 
tion (|16|l for fixed lattice spacing a — ao and for the Argon 
parameters, (c) The compressibility kt as function of tem- 
perature. 



The physical quantity (and hence k^) is obtained 
from the lowest-u intersection of the l.h.s function and of 
the r.h.s. function which is an implicit function of m^. In 
the zero temperature limit this solution corresponds to 
the harmonic limit u^iT) ~ 3kBT/ko- On the contrary, 
the second solution, which is not shown in Fig. 21 starts 
at T ~ OA' from values much higher than the nearest 
neighbors distance a and then it decreases with temper- 
ature. It corresponds to a local maximum which docs 
not minimize the free energy and it can be disregarded. 
The behaviors of Uv{T) and k'^{T) as functions of the 
temperature are plotted in Fig. [5l[a) and in Fig. [Hb). 

For a temperature independent volume, above an up- 
per temperature T*(ao) no solution is found (see Fig.|4]). 
The maximum temperature found in this way at a = ao 
for the Argon potential parameters is T* = 66.4 Uq w 
10000 K, much larger than the experimental melting tem- 
perature of the Argon T,n = 0.56 Uq ~ 82 K. Note that 
the temperature T* is so high that the linear behavior of 
Uy{T) oc T (valid for an harmonic crystal) is confined to 
a very small temperature range T ^T* (Fig. [SJa)). 

We would like to stress that, as remarked in Ref. HO, 
the breakdown of the SOFT solution at T* is not signal- 
ized by any precursor singular behavior in nor in fcv 
On this regards it is worth to analyze explicitly even the 



behavior of the compressibility kt = —{dP/dV)T/V as 
function of the temperature close to the instability tem- 
perature T* where the solution of Eq. (jl6p disappears. 
The SCPT solution for the compressibility kt is thus 
reported in Fig. ^c), showing that kt is well-behaved 
as the instability temperature at fixed volume T* is ap- 
proached. As we are going to see, the behavior of kt 
is qualitatively different when the thermal lattice expan- 
sion as a function of the temperature at fixed pressure is 
considered. 



B. Constant pressure 

Before discussing explicitly a system at constant pres- 
sure, let us consider first the interatomic distance a as 
an external and tunable parameter in order to gain some 
useful preliminary insight. 

In Fig. [5] we show the behavior of the graphical solu- 
tion of Eq. (jl6p for fixed temperature as varying a. 

We observe that the global effect of increasing the lat- 
tice spacing is to reduce the effective clastic constant fcv 
and to increase thus the mean lattice fiuctuations u^.. 
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Figure 6: Graphical solution of the self consistent equation 
Eq. (|16|l at given temperature T = 1.5 f/o and at different 
lattice spacing a > ao. The solid line and the dotted lines 
represent, respectively, the l.h.s. and the r.h.s. function of 
Eq. (|16|l for the Argon parameters. 



Thus, according to the Lindemann criterion, we can ar- 
gue that the breakdown of the soHd phase can occur at 
temperatures lower than the one corresponding to the 
constant volume case. 

In a more compelling way, the physical value of the lat- 
tice parameter a(T) (or the volume per particle V{T)/N) 
is obtained for a given pressure P and temperature T by 
inverting the relation dF^{T, V)/dV\T = -P, Eq. ([TTl). 
In Fig. [7] we show as an example the case of zero pres- 
sure P = which simply corresponds to the minimization 
of the free energy F-^ with respect to the volume V or, 
equivalently, to a. The evolution of the minimum of 
at different temperatures shows the thermal expansion of 
the system at constant pressure P = 0. 

The corresponding results for the elastic constant fcv 
and for the thermal fluctuations are reported in Fig. 
ISJa) and in Fig. ^h) . For comparison in Fig. [5] it is 
also shown the behavior of and fcv obtained by the 
self-consistent solution of Eq. p6| as functions of the tem- 
perature but assuming three representative fixed lattice 
spacing a. Note that the physical behavior of the elas- 
tic constant fcv as a function of the temperature is quite 
different in the two cases, with a decrease of fcv as a func- 
tion of T at fixed pressure whereas a hardening of fcv is 
predicted when working at fixed volume 

Our calculated data in Fig. |S] (case P = 0) extend up 
to the spinodal temperature Ts « 1.55?7oi above which 
the minimum of Pv as a function of a disappears. It 
is worth noting that such temperature Ts ~ 1.55 C/q is 
much lower than T* w 66.4 Uq related to existence of the 
self-consistent solution for the lattice fluctuation, and it 
is of the same order of magnitude of the experimental 
melting temperature Tm ~ 0.5 Uq. 

Similar results are obtained for fixed finite pressure 
P 0. As mentioned in Sec. HIl the thermal lattice ex- 
pansion can be traced by minimizing the Gibbs free en- 
ergy or, equivalently, by inverting Eq. (|17p to obtain a as 



Figure 7: Dimensionless free energy /v = F^/NUo as a func- 
tion of a in the Argon case for different temperatures. The 
filled symbols represent the minima of F^/NUo and the dot- 
dashed line shows the their evolution with T. Inset: a close- 
up for T ^Ts ^ 1.55 [/o (T = 1.35, 1.40, 1.45, 1.50, 1.55 f/o). 
For sake of visibility, in the main panel and in the inset each 
curve of the dimensional free energy has been shifted by a fac- 
tor F^/NUo -)■ F^/NUo + 3kBT/2Uo, which does not affect 
however the determination of a at the minima. 
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Figure 8: (a) mean lattice fluctuations and (b) effective 
elastic constant fcv as functions of temperature at constant 
pressure P — 0, obtained from the minimization of Pv with 
respect to a{T) as illustrated in Fig. [71 The dashed lines 
correspond to the solution of the self-consistent equations for 
fixed lattice spacing a = 1.0, 1.1, 1.2ao. 



a function of T and P. The behavior of the average fluc- 
tuation, the effective elastic constant and the lattice con- 
stant are reported in Fig. ^ Available experimental data 
of a{T) in Argon, in an experimental setup at saturated 
conditions corresponding to P w 10~^Uo/rQ « are 
also reported in FigjSJc), showing an agreement between 
our calculations and the experimental data. We observe 
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Figure 9: (a) Mean lattice fluctuations Uv, (b) effective elas- 
tic constant fcv and (c) lattice parameter a as functions of T 
for P = 0.0,0.5, 1.0[/o/?"o- The corresponding spinodal tem- 
peratures are Ts = 1.55, 1.80, 2.05 C/o- The triangles in (c) 
correspond to experimental data after Ref. [sp- 



here that increasing the pressure leads to a weaker in- 
crease of (weaker decrease of fc^) as a function of T . 
In other words, the finite pressure leads to an increase of 
the spinodal temperature Ts{P) at which the soUd phase 
disappears as a metastable solution. We also note that, 
in similar way as in the analysis at constant volume, nei- 
ther u^, fcv, nor the lattice constant a show any singular 
behavior at Ts{P)- 

In the present case at constant pressure, however, an 
additional tool of investigation is provided by the thermal 
lattice expansion coefficient a defined as: 



Vdf 



(20) 



The temperature dependence of a at fixed temperature 
is reported in Fig. I10[ showing the singular behavior at 
Ts- In the inset we plot the thermal expansion coefficient 
a as a function of the reduced variable t = \1 — T/Ts\ in 
a log-log scale, showing a power-law divergence a w 
where C ~ 0.6 was obtained by a fitting procedure. 

At the spinodal temperature, the divergence of the lat- 
tice expansion coefficient a coincides with a singular be- 
havior of the compressibility kt as it could be guessed 




kel/Uo 

Figure 10: Temperature behavior of the thermal expansion 
coefficient a/ce{T — 0) for three representative pressures P = 
0.0,0.5,1.0C/o/ro. Inset: logarithmic plot of Q/a(r = 0) vs. 
the reduced variable t — \1 — T/Ts\, showing the power law 
behavior with a ~ t~'' . 



by the general thermodynamical relationi^ 
9 C'p — C'v 



a 



VT 



-KT ; 



(21) 



where Cp and Cv are respectively the specific heat at 
constant pressure and constant volume. Note that a gen- 
eral condition for the lattice stability is the positiveness 
of the two specific heats ordered as Cp > Cy > 0.— Us- 
ing thermodynamical relations, we can rewrite Eq. (j2ip 
in the more convenient form 



a 



1 

T 



P 



dV 



(22) 



where the internal energy per particle is defined in the 
Eq. (dni), [/ = (i?)v From Eq. the divergence of 

a at the spinodal temperature Ts implies the divergence 
of Kt provided the quantity dU/dV has a finite and reg- 
ular behavior as T approaches Tg. As can be seen from 
Fig. [Til both U and its derivative dU/dV have a regu- 
lar behavior for any value of temperature and pressure, 
pointing out that a is proportional to kt at the spinodal 
temperature Ts- 



IV. PHASE DIAGRAM 

In the previous Sections, we have identified two dif- 
ferent mechanisms for the breakdown of the solid state 
phase within the SCPT model, working respectively at 
fixed volume and at fixed (moderate) pressure. In the 
first case the solid phase instability is pointed out by the 
breakdown of the self-consistent solution of the SCPT 
at a temperature T*, with no other evident signature in 
the thermodynamical properties. On the other hand, a 
solid state phase is always defined in the second case, 
which however disappears as a metastable minimum at 
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a spinodal temperature Tg, where the thermal expansion 
coefficient a and the compressibihty kt diverge. 

The competition between these two different behaviors 
can be better clarified analyzing in more details some 
typical isothermal curves in a wide range of temperatures 
as shown in Fig. [121 We can thus distinguish between 
two different pressure ranges. 

For high pressures (P > Pc with Pc « 1.02J7o/ro) we 
have only one volume solution for given pressure by in- 
creasing temperature. In this range the isothermal curves 
are monotonic as function of pressure, corresponding to a 
positive isothermal compressibility. This behavior holds 
true at high temperatures (T > 2Uo) for any pressure. 
Increasing the temperature at constant pressure leads to 
a rapid increase of the volume until the self-consistent 
solution of the SCPT disappears at the temperature T* . 
At P ~ 1.3Pc we have for instance T* = 6.50 C/q (Fig. 

ini). 

Below the critical pressure Pc, the isothermal curves 
are non-monotonic with three possible volume solutions 
for given pressure, and an intermediate region Vi{P) < 
V < V2(P) of negative compressibility kt < (dashed 
lines in Fig. [T^ . In this case, the physical more stable 
solution is the one with smaller volume4£ Increasing the 
temperature thus leads to a slight increase of V until the 
first physical solution disappears at the spinodal temper- 
ature Ts where the dP/dV = 0, and the compressibility 
Kt diverges (for P ~ 0.4, we have Ts « 1.75 Uq with cor- 
responding to a volume V « 2Vb, see Fig. [12]). At higher 
temperature, disregarding phase separated phases, the 
only physical solution is associated with the large volume 
cases (for P ~ 0.4, we have V « 8Vb for T > 1.75C/o)- 
This sudden jump of the volume at T = Ts reflects thus 
the divergence of the lattice expansion coefficient a at 
the spinodal temperature. 

Increasing further the temperature, even this high vol- 
ume solution will break due to the breakdown of the self- 
consistent SCPT solution at a high temperature T*. As 
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1 1.2 1.4 1.6 1.8 2 2.2 

V(T) / V(0) 



Figure 11: Internal energy [/ as a function of the volume 
y(r) varying the temperature for different pressures P = 
0.0, 0.5, 1.0(7o/rg. 



a matter of facts, this kind of instability is driven by the 
disappearing of the SCPT solution at T* (the same en- 
countered when working at constant volume) in contrast 
to the compressibility divergence at Ts- 

We can summarize all the above discussion in a com- 
pact phase diagram. Fig. I13[ which reproduces on a 
quantitative level the schematic phase diagram sketched 
in Fig. [2] 

The different kinds of instabilities of the solid phase as 
a function of the pressure are depicted in Fig. 1141 where 
the dashed line represents the spinodal temperature Ts 
accompanied by the divergence of the lattice expansion 
coefficient a and of the compressibility kt , while the solid 
line marks the disappearing of the SCPT solution, as 
function of temperature, not reflected in any singular be- 
havior of the thermodynamical properties. Note that for 
T < Ts, T* is not defined. 

The physical mechanism of the breakdown of the solid 
phase is determined by the first instability encountered 
by increasing the temperature. For P < Pc this is given 
by the compressibility divergence at Ts- For P > Pc 
no thermodynamical instability is encountered and the 
unique instability of the solid phase is associated with 

Fig ll4l permits also to clarify in simple terms why the 
mechanical instability at constant volume occurs at much 
higher temperature T* and with different phenomenology 
than at constant pressure. Considering for simplicity a 
constant volume a = qq {V — Vq), increasing the temper- 
ature is indeed associated with an effective increase of P. 
Such increase is however negligible on the scale P « Pc, 
so that the first encountered solid phase instability oc- 
curs at much higher temperatures, as shown in the inset, 
and it is related to the SCPT breakdown at T* . 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

P/Pc 

Figure 12: (color online) Representative isothermal curves of 
different regimes, for T = 0, 1.55, 1.75, 2.05, 3.55, 5.00, 6.50 Uo- 
Solid line regions represent the case where the compressibility 
Kt is positive, while dashed line regions correspond to kt < 
0. The volume is expressed in units of the Vb = Nuq, and 
the pressure in units of the critical pressure Pc = 1.02 Uo/rQ, 
whose meaning is discussed in the text. 
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12 3 4 



Figure 13: (color online) Phase diagram for the solid insta- 
bilities in the volume-temperature space. The SCPT solution 
disappears for T > T* . The long dashed line represents the 
spinodal temperature Tg where the compressibility is nega- 
tive (kt < blue shadow region) for T < Ts- Also shown 
is the isobaric curve (tiny solid black line) corresponding to 

P = Pc. 




P / Pc 

Figure 14: Pressure dependence of Ts and T* for the Argon 
case. Also shown as symbols is the P vs. T dependence 
of the isochoric line for V = Vo- Inset: same quantities in a 
large axes scale, in order to show the intercept of the isochoric 
line with T* which defined the solid phase instability at fixed 
volume V = Vq. 



V. DISCUSSION AND CONCLUSIONS 

In this paper wc have revisited the self-consistent 
phonon theory (SCPT) to analyze the mechanism of the 
mechanical instability of an overheated metastable solid 
at finite pressure. 

Within this approach, we have predicted that two dif- 
ferent instabilities are possible for the mechanical melt- 
ing. They occur at very different temperatures, depend- 
ing whether working at fixed volume or pressure. In the 
first case the instability of the solid phase occurs at rather 
high temperatures and it is associated with a breakdown 
of the SCPT solution at T*, not reflected in any sin- 
gular behavior of the thcrmodynamical quantities. On 



the other hand, the mechanical instability approached 
at fixed pressure, for pressures smaller than a critical 
pressure Pc, can be properly interpreted in terms of a 
spinodal temperature Ts at which the solid phase, al- 
ways defined, disappears as a metastable minimum. In 
the latter case the mechanical instability is refiected in 
the divergence of the thermal compressibility kt as well 
as in the divergence of the lattice expansion coefficient a 
at r = T5. Such kind of spinodal instability disappears 
for P > Pc, where only the instability at T* driven by 
the breakdown of the SCPT solution remains. The value 
of Pc is found to be of order Pc ~ Uo/aQ where Uq is the 
energy minimum of the pair potential and oq the lattice 
spacing. For the weakly bound rare-gas solids we have Uq 
of the order of few meVs, corresponding thus to a critical 
Pc of the order of some MPa. 

For a quantitative comparison with real systems, in 
this paper we have mainly focused on the specific case of 
an Argon classical solid where intensive theoretical and 
experimental investigation have been performed in liter- 
ature (as example, see Ref. [2^ and references therein). 

Similar analysis can be however carried out in a simple 
way as well for other rare-gas solids (Xe, Ne, Kr). The 
results are collected in Table |T] where we report the the- 
oretical T* and Ts compared also with the experimental 
melting temperature T,„. As we can see, Tg is system- 
atically much lower than T* and it is of the same order 
of magnitude of the experimental melting temperature 
T„j. Inclusion of higher order anharmonic terms^ and 
the development of models beyond the Einstein one^i^ 
(i.e. taking into account the full phonons dispersion) 
might further reduce Tg towards the empirical range for 
overheating T^ax 1.5T„3 





Tm (K) 


T* (K) 


Ts (K) 


Ne 


25 


2450 


69 


Ar 


83 


9747 


228 


Kr 


116 


18986 


317 


Xe 


161 


21662 


445 



Table I: Experimental melting temperature Tm of (Ne, Ar, Kr, 
Xe) compared to the temperatures of the mechanical lattice 
instability as obtained by the SCPT at fixed volume (T*) with 
a — ao and at fixed pressure (Ts) with P = . 

Such analysis points out that taking into account the 
lattice expansion is of fundamental importance not only 
for a proper quantitative estimate of the instability tem- 
perature of the solid phase, but also for revealing the 
different instability mechanisms associated or not asso- 
ciated with a singular behavior of the thcrmodynamical 
properties. 
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